Quenching of high-p T hadrons: Energy Loss vs Color Transparency 
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High-pT hadrons produced in hard collisions and detected inclusively bear peculiar features: (i) 
they originate from jets whose initial virtuality and energy are of the same order; (ii) such jets 
are rare and have a very biased energy sharing among the particles, namely, the detected hadron 
carries the main fraction of the jet energy. The former feature leads to an extremely intensive gluon 
radiation and energy dissipation at the early stage of hadronization, either in vacuum or in a medium. 
As a result, a leading hadron must be produced on a short length scale. Evaluation within a model 
of perturbative fragmentation confirms the shortness of the production length. This result is at 
variance with the unjustified assumption of long production length, made within the popular energy 
loss scenario. Thus we conclude that the main reason of suppression of high-pr hadrons in heavy ion 
collisions is the controlled by color transparency attenuation of a high-pr dipole propagating through 
the hot medium. Adjusting a single parameter, the transport coefficient, we describe quite well the 
data from LHC and RHIC for the suppression factor Raa as function of pr, collision energy and 
centrality. We observe that the complementary effect of initial state interaction causes a flattening 
and even fall of Raa at large pr- The azimuthal anisotropy of hadron production, calculated with 
no further adjustment, also agrees well with data at different energies and centralities. 



I. INTRODUCTION 

A colored parton produced with a high momentum in 
a hard reaction hadronizes, forming a jet of hadrons. It is 
natural to expect that the production time of such a jet, 
averaged over jet configurations, rises with the jet energy 
due to the effect of Lorentz time dilation. Although the 
jet is detected at macroscopic distances from the collision 
point, its space-time development at the early stages of 
hadronization can be probed with nuclear targets 

In this paper we concentrate on the rare type of jets 
in which the main fraction of the jet momentum is car- 
ried by a single (leading) hadron. In some cases, like 
in semi- inclusive deep-inelastic scattering (SIDIS), such 
events can be selected explicitly, because the fractional 
light-cone momentum of the detected hadron can be 
measured. In high-p^ single hadron production processes 
the fractional hadron momentum is unobserved, but the 
convolution of the steeply falling jet momentum distribu- 
tion with the fragmentation function picks up high values 
of Zh (see Sect. Ill Dj) . Thus, inclusive high-p^ hadron pro- 
duction without observation of the whole jet implicitly 
selects an unusual type of jets with a very biased sharing 
of energy. 

Another generic feature of such jets is an extremely 
high initial virtuality, which is of the same order as the 
jet energy. This leads to a very intensive gluon radiation 
and energy dissipation at the early stage of hadroniza- 
tion. In order to respect energy conservation in the pro- 
duction of a high-z/j hadron, the radiative dissipation of 
energy must be stopped by the production of a colorless 



hadronic configuration (QCD dipole or pre-hadron) on a 
short time or length (we use both terms interchangeably) 
scale. This is considered in detail in Sect. UH where the 
rate of radiative energy loss in vacuum is calculated as 
a function of time. The production length l p of a col- 
orless dipole finalizing hadronization is calculated within 
a model of perturbative hadronization and found to be 
rather short. The important observation is a weak de- 
pendence of l p on pr, which might look counterintuitive, 
because the Lorentz factor is expected to stretch l p at 
higher px- However, the rate of energy dissipation in- 
creases as well, and this works in the opposite direction, 
trying to shorten l p . 

Since a colorless dipole is created at a short time scale 
inside a dense medium, it has to survive through the 
medium in order to be detected. The evolution of the 
dipole in the medium and its attenuation is the subject 
of Sect. IIIII The key phenomenon controlling the dipole 
surviving probability is color transparency, which cor- 
responds to the enhanced transparency of the medium 
for small-size dipoles 0. We employ the relation be- 
tween the dipole cross section and transport coefficient 
(broadening) found in 3-5]. Correspondingly, the ob- 
served magnitude of hadron attenuation can be used as 
a probe for the transport coefficient, which characterizes 
the medium density. 

In Sect. lIVI we compare the calculated suppression fac- 
tor Raa with data, as function of pt, collision energy and 
centrality of the collision. This comparison involves only 
one fitted parameter, g . which is the maximal transport 
coefficient of the medium created in a central collision of 
given nuclei, at a given energy. Otherwise, this param- 
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eter is universal for all observables. The shape of the 
Pt dependence of Raa is found to be in a good accord 
with data. In particular, the observed rise of Raa(.Pt) at 
LHC is easily and naturally explained by the color trans- 
parency effect, calculated within the rigorous quantum- 
mechanical description known as the path-integral tech- 
nique. 

Comparison with data results in the transport coef- 
ficient , which ranges from q = 1.2GeV 2 /fm at ^/i = 
62GeV up to 2GeV 2 /fin at y/s = 2.76 TeV, for colli- 
sions of heavy nuclei, gold and lead. These values of the 
transport coefficient are about twice as large as those 
that were found in Q within a simplified model of dipole 
evolution. Nevertheless, they are an order of magnitude 
smaller than what was found in the analysis Q , based on 
the energy loss scenario (see e.g. @), which relies on the 
unjustified assumption of a long production length l p . 

It is worth emphasizing that our approach, based on 
perturbative QCD, is irrelevant to data at pr ^5 6GeV, 
which are apparently dominated by hydrodynamics. 

An additional effect related to initial state interactions 
(ISI) of the colliding nuclei is described in Sect. IIVBI 
The excitation of higher Fock states by multiple inter- 
actions leads to enhanced nuclear suppression of particle 
production with large xt and/or xl [€>}■ This effect can 
be seen in the px dependence of Raa at the RHIC en- 
ergies Vs = 200 GeV and 62 GeV. Also LHC data at 
yfs = 2.76 TeV indicate a leveling of the Raa behavior 
at the maximal measured px, and we expect even a fall 
at p T > 100 GeV. 

A complementary test of our approach is provided by 
data on azimuthal anisotropy of produced hadrons, as is 
described in Sect. IIV CI In fact, we explain the measured 
difference between Raa for in- and out-of-plane events 
and the asymmetry parameter V2(pt), with no additional 
adjustment. 



II. ENERGY CONSERVATION: HOW LONG 
DOES HADRONIZATION LAST ? 

One should discriminate between the observation of a 
jet initiated by a parton produced in a hard reaction (e.g. 
high-pT processes or deep-inelastic scattering (DIS)), and 
the detection of only a single hadron produced inclusively 
with a large fractional light-cone momentum in a hard 
process at high energies. The latter process corresponds 
to a very rare jet configuration, where the main frac- 
tion Zh of the jet energy E is carried by a single hadron, 
while all other hadrons in the jet must share the smaller 
energy (1 — Zh)E. The deficit of energy imposes certain 
constraints on the space-time development of such a jet, 
which is different from an averaged jet. Our definition of 
the characteristic time scales is illustrated in Fig. [TJ The 
quark regenerating its color field, which has been stripped 
off in a hard reaction, intensively radiates gluons and dis- 
sipates energy, either in vacuum or in a medium. Multi- 
ple interactions in the medium induce additional, usually 
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FIG. 1: (Color online) Space-time development of hadroniza- 
tion of a highly virtual quark producing a leading hadron, 
which carries the main fraction Zh of the initial quark light- 
cone momentum. 



less intensive, radiation. The loss of energy ceases at the 
moment, call production time t p , when the quark picks 
up an antiquark neutralizing its color. The produced 
colorless dipole (also called pre-hadron) does not have 
either the wave function or mass of the hadron, but it 
takes the formation time tf to develop both. The forma- 
tion stage is described within the path-integral method 
in Sect. ITTlBl 



A. String model 

The simple example of the consequences of energy con- 
servation is the string (or color flux tube) model [10( . The 
usual expectation is that the higher is the energy of the 
produced hadron, the longer it takes to be produced, as 
follows from Lorentz time dilation. On the contrary, in 
11] it was found that the production time of a hadron 
with energy z^E shrinks down to zero upon reaching 
the kinematic bound of maximal energy E. Indeed, the 
quark that initiated a jet is losing its energy with the 
rate dE/dt = —re, where k w 1 GeV/ fm is the string 
tension. The energy loss comes from the hadronization, 
which is developed as a series of string breaks by qq pairs 
tunneling from the vacuum (Toj . 

Since the leading quark keeps losing energy, it can pro- 
duce a hadron with energy ZhE only within a certain time 
interval, restricted by energy conservation, 
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(1) 



Such a shrinkage of the production time was explicitly 
confirmed in Montecarlo models [HI, [H| . 

The string model is mentioned here only as an example 
of the constraints imposed on t p by energy conservation. 
It should not be taken literally as a hadronization mech- 
anism for hard reactions, where a highly virtual parton 
neutralizes its color perturbatively on a short time scale 

B 

At this point we should specify our terminology in or- 
der to avoid further confusions. Indeed, different time 
scales are debated in the literature, coherence time, for- 
mation time, production time. What we call produc- 
tion time (in Eq. ([1} and in what follows), is the time of 
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color neutralization of the leading quark by an antiquark 
picked up from the string or generated perturbatively 
[14| . Notice that this is not yet the final hadron, which is 
characterized by both a specific wave function and mass. 
What is produced at the time scale t p is a colorless qq 
dipole, having certain separation, but not mass, which 
we call conventionally pre-hadron. It takes the forma- 
tion time to develop the wave function, 



2E h 
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(2) 



where Eh — z^E; rrih and m/j* are the masses of the 
hadron (we assume it to be the ground state) and the first 
radial excitation h*. This time scale does not shrink at 
Zh — > 1, but keeps rising. It can be derived in the multi- 
channel representation for the diffractive scattering as 
the inverse minimal value of the longitudinal momentum 
transfer in the off-diagonal diffractive transitions between 
different states. 

Eq. @ can also be understood in terms of the un- 
certainty principle. Namely, the produced qq dipole has 
a certain size and can be projected either to h, or to 
h* . According to the uncertainty principle, it takes the 
proper time tj — l/(m/j» — m/j) to resolve between these 
two levels. Applying the Lorentz boosting factor we get 
©■ 

Concerning the coherence time scale, this is a more 
general term, which means that quantum-mechanical in- 
terferences are important. Depending on the context and 
the employed theoretical tools, it might play role on ei- 
ther the production time or the formation time. 



B. Radiational energy loss 

As a result of hard interaction characterized with the 
scale Q 2 , the parton is produced with part of its color 
field stripped off, up to transverse frequencies kx < Q- 
Hadronization of such a highly virtual quark cannot be 
described adequately in terms of the nonperturbative 
string model. Regeneration of the quark color field is 
associated with radiation of gluons, which take away a 
part of the quark energy and contribute to the forma- 
tion of the jet. In fact, at high virtualities Q 2 this gluon 
radiation becomes the dominant source of energy loss in 
vacuum. 

One should strictly discriminate between vacuum and 
medium-induced energy loss. The former includes the 
lost energy, which goes into gluon radiation and/or into 
setting up the string field, in other words into jet for- 
mation. The latter corresponds to the additional energy 
loss caused by the multiple interactions of the jet in the 
medium. The vacuum rate of energy loss usually signif- 
icantly exceeds the medium-induced one. Here we con- 
centrate on the study of the hadronization pattern in 
vacuum. 

The time dependent radiational energy loss can be cal- 



culated as follows [l], [H, HH 

Q 2 l 
&E rad (t) = E J dk 2 J dx : 
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dx dk'- 



■ e(t-t»), (3) 



where the coherence time for radiation of a gluon with 
fractional light-cone momentum x and transverse mo- 
mentum k reads, 



t'i 



2Ex(l 



k 2 + x 2 m 2 



(4) 



The step function in Eq. (j3|) excludes from the integra- 
tion those gluons which are still in coherence with the 
radiation source, and did not materialize on mass shell 
during the time interval t. The soft cutoff A in Eq. ([3]) 
is fixed at A = 0.7 GeV. This choice is dictated by data 
(see in [22|, [23| ) , which indicates a rather large primordial 
transverse momentum of gluons. 

The spectrum of radiated gluons in Eq. ([3]) has the 
form, 



dn„ 



2a s (k 2 ) fc 2 [I + (I-a;) 2 
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(5) 



where a s (k 2 ) is the running QCD coupling. 

The time dependence of radiational energy loss in vac- 
uum exposes a nontrivial behavior [l|, [T^ . During the 
time interval t < (2E / Q 2 )(l — Zh) the energy loss rises 
linearly with time, 



AE rad (t)=t^(Q 2 



A 



(6) 



However, at larger t the rate of energy loss starts falling, 
AE rad (t) is leveling off and at t > (2E/A 2 )(1 - z h ) the 
gluon radiation completely ceases. Then the quark loses 
energy only via nonperturbative mechanisms (strings). 

Similar to Eq. (JTJ) , the production time is restricted by 
energy conservation to, 



AE{t p ) « E(l - z h ). 



(7) 



Apparently, the increase of the energy loss rate in 
Eq. ([TJ, k K + \dE rad /dt\, caused by gluon radiation, 
can only shorten the production time. 



C. Peculiar aspects of high-fc^ jets 

The solution of Eq. ([7]), the production time of a lead- 
ing (pre)hadron, depends on the jet energy and virtuality. 
In deep-inelastic scattering (DIS) these are two indepen- 
dent variables, and usually E 2 3> Q 2 . 

For a parton produced with high transverse momen- 
tum kr normal to the collision axis in the cm. frame, 
its energy E kx- The hard scale for such a process is 
also imposed by the transverse momentum, i.e. Q 2 = k T . 
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Thus, a high-Zcr jet is in an unique kinematic domain of 
extremely high virtuality, Q 2 = E 2 , which cannot be ac- 
cessed in DIS. This fact leads to a specific behavior of the 
production time for high-p-r hadrons, which is different 
from what is usually measured in SIDIS. 

Indeed, we can trace the dependence of t p on E and 
Q 2 using the approximate relation analogous to Eq. ([TJ), 



t < 



E 



(\dE/dt\) 



(l-z h ). 



(8) 



If one increases the jet energy keeping the virtuality Q 2 
fixed, so that (\dE/dt\) does not vary, the production 
time rises linearly with E and eventually exceeds the 
time of jet propagation through the medium. Indeed, 
the observed magnitude of nuclear suppression of lead- 
ing hadrons in SIDIS decreases with energy and vanishes 
at jet energies E ~ 100 GeV (2(j. However, at medium 
and low energies [H, HU the energy loss scenario [l7j . 
which assumes a long production time, fails to describe 
the nuclear suppression observed in SIDIS at large z%. 

And vice versa, if one keeps the energy E fixed, but 
increases the virtuality Q 2 , the mean rate of energy loss 
in the denominator of (|SJ rises, and the production time 
shrinks. 

Therefore, it is not obvious what will happen to t p if 
both the energy and virtuality rise simultaneously, as it 
happens for high-fcr jets. In spite of the rising Lorentz 
factor, the jet virtuality Q 2 = E 2 rises as well and causes 
a dramatic enhancement of radiative energy loss, which 
may result in a shorter production time. 

Eq. describes the time dependence of the energy 
radiated by a virtual quark. The total amount of radia- 
tional energy loss is given by the same equation without 
the 9-function, 



AE t , 



3vr \ A 



(9) 



where we assumed x -C 1 and fixed a s for the sake of 
simplicity. Let us inverse the problem of the time depen- 
dence of energy loss Eq. ((3]) and ask: how long does it 
take for a highly virtual quark or gluon with Q 2 = E 2 to 



radiate a fraction 



6(t) 



AE(t) 
AE tot 



(10) 



of the total radiated energy? Solving Eq. one gets 



5|ln(f) if*. 1 I* v 



(ii) 



From the second line of this equation we conclude that 
a high-/cT parton radiates half of the total energy loss 
during a very short time interval, t(S = 1/2) = 2/(eA) 
0.2 fm. This interval is independent of energy. 
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FIG. 2: The path length taken by a parton (quark or gluon) 
of energy E and virtuality Q 2 = E 2 to radiate a fraction 
S — AE/E of the total vacuum energy loss. The curves cor- 
respond to different jet energies, E — 10, 20, 50 and 100 GeV. 



The ^-dependence of the time interval t(S) is illustrated 
in Fig. [5] for several jet energies. This confirms that 
a high-fc^ parton radiates the main fraction of the en- 
ergy vacuum loss at the early stages of hadronization, on 
very short time intervals. Notice that inclusion of non- 
perturbative mechanisms of energy loss should lead to 
even faster dissipation of energy. Thus, we conclude that 
the fast degradation of energy of a highly virtual parton 
makes impossible the production of leading (pre)hadrons 
on a long time scale. 

So far we explored the time dependence of radiation 
with no constraints on the radiation process. However, 
detection of a hadron carrying a large fraction z% of the 
initial parton energy essentially affects the energy loss 
pattern due to energy conservation. The production of 
leading hadrons in such jets was studies in [21j , where it 
was found to have the form, 



dt v 



oc (1 - z h )S(t p ,z h ) 



(12) 



Here Q 2 = k 2 ^, and the time dependent fractional mo- 
mentum Zh(t) reads, 



Zh{tp) 



Zh\ 
X I 



Zh 1 + 



AE(t P/ 
E 



+ [z h {l-z h f 



(13) 

The energy loss AE(t p ) includes both perturbative and 
nonperturbative (strings) sources of energy dissipation. 
The former is calculated with Eq. in which an ad- 
ditional kinematical constraint is introduced: the energy 
of radiated gluons, cj = aE + k 2 /4aE, should not ex- 
ceed the bound uj < (1 — Zh)E. Such a ban for radiation 
of part of the gluon spectrum during the time interval 
t < t p leads to a suppression known as Sudakov fac- 
tor S(t p ,Zh), introduced in Eq. (fT2|). It is defined as 
S(t p ,Zh) = exp [— (n g (t p , Zh)}], where the mean number 
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of gluons (n g (t p ,Zh)} which have radiation time shorter 
than t p , but cannot be radiated due to energy conserva- 
tion, is calculated with the same gluon spectrum Eq. (|5|) . 
Examples for the Sudakov factor at different values of 
Zh are shown in Fig. [3j as functions of the production 
length lp = t p and for jet energies E = 10 and 100 GeV. 
We see that energy conservation vetoing part of the radi- 
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FIG. 3: Sudakov suppression factor caused by the ban for 
radiation of gluons with fractional energy higher than 1 — zh, 
during the time interval t < t p = l p . Calculations are done 
for jet energies E — Q — kr — 10, 100 GeV, and several 
fractional hadron momenta z^ = 0.6 — 0.9. 

ation spectrum results in the Sudakov suppression factor, 
which substantially reduces the production time Eq. (|12|l . 

Now we are in a position to perform numerical cal- 
culations for the production time distribution func- 
tion. Examples for the differential fragmentation func- 
tion Eq. ([12]). at Q 2 = E 2 = v\l z \i are depicted in 
Fig. 0]as function of l p . Since the absolute value of the 
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FIG. 4: The differential fragmentation function Eq. (|12p (in 
arbitrary units) at Q 2 = E 2 as function of l p for quark jets 
with energies E = kr — 10, 20 and 100 GeV (from top to 
bottom) and Zh = 0.7. 



fragmentation function steeply varies with Zh, we renor- 
malized it by adjusting to the same value at l p = 0, and 
plot it in arbitrary units. 

As a test of the modeled differential fragmentation 
function Eq. (I12p , we integrated it over t p and compared 
with data. Our result reproduces pretty well [2l| the phe- 
nomenological function Dt T / q (zi 1 ,Q 2 ) [25| fitted to data, 
at Zh ^ 0.5. Correspondingly, we consider our calcu- 
lated tp-distribution Eq. (fl"2"j) to be trustable within this 
interval of z^. 

Eventually, using the distribution (fT2j) we can calculate 
the mean production time, 



/* < 1 f M t dD v/q (z h ,E 2 ) 



dt r 



(14) 

The results are presented in Fig. [S] Naturally, the pro- 
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FIG. 5: (Color online) The mean production length as 
function of energy for quark (solid curves) and gluon(dashed 
curves) jets. In both cases the curves are calculated at 
Zh = 0.5, 0.7, 0.9 (from top to bottom). 

duction length for leading hadrons in jets initiated by 
gluons is shorter than for quarks, because the dissipation 
of energy in gluon jets is more intensive. 

Notice that color neutralization, resulting in the pro- 
duction of a pre-hadron, also may be subject to coher- 
ence. This means that the amplitudes with different t p 
can interfere, so one cannot identify with certainty the 
production moment t p . One may be even unable to say 
whether the pre-hadron was created inside or outside the 
medium. Such a quantum-mechanical uncertainty was 
explicitly demonstrated in a SIDIS example [32j . where 
the interference term in the production cross section was 
found to be 100% important and negative. However, in 
the case of a high-py jet the starting virtuality is so high, 
and the dissipation of energy so intensive, that every am- 
plitude is constrained to have a short t p . Therefore, in 
what follows we rely on the probabilistic description of 
the space-time development illustrated in Fig. [1] 
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D. The mean value of Zh 

As expected, the production time varies with Zh, which 
unfortunately cannot be measured in the process under 
consideration, but one can evaluate its mean value and 
then rely on it in further calculations. As was stressed 
above, inclusive production of hadrons with large trans- 
verse momentum enhances the large-z^ part of the frag- 
mentation function D(zh,Q 2 )- This happens due to the 
steepness of the fey-spectrum of the produced partons, 
quarks or gluons, which has to be convoluted with the 
fragmentation function. This convolution defines the 
mean fractional momentum (z/J. 

First of all, we should check how well we can describe 
data in pp collisions. We employ the simple model pro- 
posed in [331 ] . based on fey-factorization. 



- — = K / dxidxjd 2 k i T<i 2 k i T 

''!<•' i" ax, J 

x Fi/ p (xi,kiT, Q 2 ) Fj/ p (xj,kj T , Q 2 ) 

x ^kl) — D h/k {z h ,Q 2 ). (15) 

dt 7T Zh 

Here da(ij — > kl)/dt is the cross section of parton scat- 
tering; the kinematic variables and their relations can be 
found in [33[ . Following [33| we assume a factorized form 
of the transverse momentum distribution, 
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F i/p (x,k T ,Q 2 ) = F i/p (x,Q 2 )g p (k T ,Q 2 ), 



where 



g p {k T ,Q 2 ) 



tt(4(Q 2 )) 



-4/(4(q 2 )) 



(16) 



(17) 



The scale dependence of (fe^(Q 2 )) was parametrized in 
33] as (k 2 ) N (Q 2 ) = 1.2 GeV 2 + 0.2a s {Q 2 )Q 2 , with pa- 
rameters adjusted to next-to-leading order calculations. 

We use the phenomenological parton distribution func- 
tions (PDF) F i/p (x,Q 2 ) from MSTW08 leading order 
(LO) [34|. For the fragmentation function D h / k (zh, Q 2 ) 
we rely on the LO parametrization given in [35[ . 

The results of the differential invariant cross section 
calculations Eq. (I15[) arc compared with data at ^fs = 
200 GeV and 7TeV [13 in Fig. ® We see that the 
employed model reproduces quite well the shape of the 
measured cross section up to the maximal available mo- 
menta. The absolute normalization (not important for 
us) is regularized by the K-factor in (|15[) . which we found 
to be k « 1 — 1.5 depending on the energy. 

Notice that while at the LHC energy the cross section 
is fully dominated by gluon jets, at y/s — 200 GeV quarks 
are important and even dominate towards the upper end 
of the available range of pr- 

Eventually, we can average the fractional hadron mo- 
mentum Zh weighted with the convolution Eq. (|15j). The 
results are depicted in Fig. separately for quark and 
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FIG. 6: (Color online) p-r-dependence of pion production 
in pp collisions at ^/s = 200 GeV (upper panel) and 7TeV 
(lower panel). The contributions of quark and gluon jets are 
shown in the upper panel by dashed and dotted curves re- 
spectively. Data points are from the PHENIX [H and CMS 
[37j experiments. 
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FIG. 7: (Color online) The mean fraction (zh) of the jet en- 
ergy carried by a hadron detected with transverse momentum 
Pt- The calculations are performed as described in the text, 
for collision energies = 200, 2760 and 7000 GeV. 
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gluon jets (upper and bottom solid curves) and at dif- 
ferent energies, yfs = 200, 2760 and 7000 GeV. We see 
that the lower is the collision energy, the larger is (zh), 
especially at large pr, because the parton kx distribution 
gets steeper. In the energy range of the LHC the magni- 
tude of (zh) practically saturates as function of y/s and 
Pr- 

Such a large value of the fractional jet energy carried 
by a single inclusively detected high.-px hadron (without 
observation of the whole jet), makes its space-time de- 
velopment very different from an usual high-p^ jet, when 
the whole jet is reconstructed. In the latter case, if no 
special selection is made, the mean fractional momenta 
of hadrons are very small. Correspondingly, energy con- 
servation does not impose any severe constraints on the 
hadronization time scale, which rises with px and may 
be long. 



III. ATTENUATION OF LEADING HADRONS 
IN A DENSE MEDIUM 

In previous sections we found that in rare events, in 
which the produced hadron carries the main fraction Zh 
of the initial parton light-cone momentum, the intensive 
gluon radiation and energy dissipation in vacuum by a 
highly virtual parton produced with high kx does not 
leave much time for the hadronization process. Particu- 
larly, such rare events are selected by detecting a hadron 
with large px- Fig. [7] shows that the detected hadron 
carries on average more than half of the jet momentum. 
In order to respect energy conservation, the intensive 
dissipation of energy by the parton should be stopped 
promptly by color neutralization, i.e. production of a 
colorless pre-hadron, otherwise the leading parton, which 
lost too much energy, will be unable to produce a hadron 
with large Zh- Other possibilities for reduction of the vac- 
uum energy loss by non-radiation of gluons are strongly 
Sudakov suppressed. As a result, the time scale for pro- 
duction of a colorless dipole is rather short and does not 
rise with px, as is depicted in Fig. [5] 

If this process occurs not in vacuum, but in a dense 
medium, multiple interactions of the parton generate 
more energy loss, which makes the production time even 
shorter. The further interactions of the dipole in the 
medium are mostly inelastic color exchange collisions. 
Indeed, the cross section of inelastic interactions is pro- 
portional to the dipole separation squared, r T , while the 
elastic scattering cross section is oc r T , which is negligibly 
small. So, any inelastic interaction of the colorless dipole 
with color exchange will resume the gluon radiation and 
dissipation of energy. Of course color neutralization may 
happen again via creation of a new dipole, but such a 
reincarnation of the pre-hadron will result in a substan- 
tial reduction of its momentum. 

Thus, we should evaluate the survival probability W, 
i.e. the chance for a dipole to escape from the dense 



medium having no inelastic interaction on the way out. 
Apparently, this is subject to the effect of color trans- 
parency 0], i.e. the rate of attenuation of small-size 
dipoles vanishes quadratically with the dipole transverse 
separation r, 



dW 



= -e(l)r 2 



(18) 



r T ->-0 



where e(l) characterizes the medium, and varies with co- 
ordinates and time. 

To avoid terminological confusions, notice that in high- 
Pt hadron production the momentum px is meant to 
be transverse to the collision axis. However, the "trans- 
verse" dipole separation is meant to be transverse relative 
to the vector px- 

It was found in Q that a parton propagating through 
a nuclear medium experiences broadening, whose magni- 
tude is controlled by the small-r behavior of the dipole- 
nucleon cross section. Thus, the broadening and attenu- 
ation rates in a medium turn out to be related as, 



(19) 



Here the broadening rate q(l) — dAq 2 /dl is usually called 
transport coefficient and used as a characteristics of the 
medium [27| ■ Notice that q is the transverse momen- 
tum of the parton relative to its initial direction. The 
transport coefficient is proportional to the medium den- 
sity, which is function of impact parameter and time. In 
what follows we rely on the popular, although poorly jus- 
tified model for q, which is assumed to be proportional 
to the number of participants and gets diluted with time 
as p(t) — l/t, due to the longitudinal expansion of the 
produced medium. Correspondingly, the transport co- 
efficient depends on impact parameter and time (path 
length I = t) as [28| . 



q(l,b,r) 



tpart 



(0,0) 



e(j-j<>), 



(20) 



where b is the impact parameter of nuclear collision, f is 
the impact parameter of the hard parton-parton collision 
relative to the center of one of the nuclei,, n part (b, f) is the 
number of participants, and go is the rate of broadening 
of a quark propagating in the maximal medium density 
produced at impact parameter r = in central collisions 
(b = 0) at the time t — to — h after the collision. The 
corresponding transport coefficient for gluons should be 
9/4 bigger. The equilibration time to is model dependent. 
Our results are not very sensitive to it, and we fix it at 
to = = 0.5 fm. 

As far as the process of high-p^ hadron production is 
considered as a probe for the medium properties, we treat 
the transport coefficient q(l,b,f) as an adjustable quan- 
tity. Once the shape of this function is fixed by the model 
Eq. ([2H|) . the only fitted parameter is the maximal value 
of the transport coefficient qo, which is A-dependent. 
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A. Evolution and attenuation of a dipole: 
heuristic description 

Here we employ a simplified description of the time 
evolution, in terms of the mean dipole transverse separa- 
tion. The dipole produced with a very small initial size 
r ~ starts expanding with a speed given by the 

uncertainty relation dr/dtoc 1/r [3, [29|, [38|] . Correspond- 
ingly, the /-dependence of r is described by the following 
linear differential equation , 



1 



dr 



dl r(l)Ehu(l — a) 



(21) 



where a and 1 — a are the fractions of the dipole light- 
cone momentum carried by the quark and antiquark, and 
Eh = Pt is the energy of the dipole. 
The solution of Eq. (|2"T1) reads, 



21 



a(l — a)pr 



(22) 



where tq is the initial dipole size. If it is small, r$ ~ 1/fc-r, 
its magnitude is quickly "forgotten". Indeed, the first 
term in (|2"21 starts dominating at / 3> 1/pr, and the 
value of the initial size ro is unimportant. 

The mean value of r 2 {l) Eq. (f2"2"j) can be used in 
Eq. (|T5| to evaluate the attenuation of a dipole, initially 
small and evolving its size in a medium characterized 
with a transport coefficient q. 



RAB(b, T,pr) 



27 6XP 



-— ldllq(l,b,T + j) 
p T 

l 



This is the medium attenuation factor for a dipole pro- 
duced in a hard internal NN collision at impact parame- 
ter r , in a collision of nuclei A and B with impact param- 
eter b. The produced hadron is detected at azimuthal an- 
gle <j> relative to 6, i.e. l-b = lb cos cj). The bottom limit of 
/-integration is L = max{/ p , Iq}. The dependence of the 
transport coefficient on coordinates is given by Eq. (l20t . 
In (|2"3")l we fixed a — 1/2, because the dipole distribution 
amplitude over r and a is projected to the hadron (pion) 
wave function, which has a maximum at a = 1/2 [3^ |. 
Moreover, we rely on Berger's approximation [ioj], which 
fixes a at this value. 

This simplified approach was employed in Q and de- 
scribed quite well the first data from the ALICE ex- 
periment [4l| for central collisions. Important observa- 
tions made in Q are: (i) the value of the parameter go 
needed to explain the observed suppression is an order 
of magnitude smaller than the results of the data anal- 
ysis [8J, based on the energy- loss scenario. At the same 
time, it agrees well with the perturbative evaluation of 



Qo [221 j and with the analysis [23, |38| of data on suppres- 
sion of J/ty produced in central gold-gold collisions at 
\/i = 200 GeV. (ii) The observed rising py-dependence 
of Raa{pt) [HI is naturally explained by the effect of 
color transparency Eq. (fT8|) . Namely, the higher is the 
dipole energy E^ = px (in the medium rest frame), the 
more is the Lorentz dilation of the dipole size expansion, 
and therefore the less is the absorption. This is why there 
is a factor 1/pr in the exponent of Eq. ([23)) . 

As was commented above, this simplified heuristic 
model allows to easily understand the main features of 
the underlying dynamics. It also allows to speed up the 
calculations. However, it misses some details which may 
affect the result. In particular, Eq. (l2~Tj) describes the ex- 
pansion of the dipole in vacuum. However, in a medium, 
color filtering effects modify the path-length dependence 
of the mean dipole separation. Namely, dipoles of large 
size are strongly absorbed, while small dipoles attenuate 
less. Correspondingly, the mean separation in a dipole 
propagating in a medium should be smaller than that 
given by the differential equation (|2~T|) for the dipole ex- 
pansion in vacuum. Introducing an absorptive term we 
arrive at a modified evolution equation, 



dP 
~~dT 



r 4 (Z)e(Z), 



(24) 



E h a(l - a) 

where e(l) is related by (fl~9|) to the transport coefficient. 
Apparently, such a modification results in a reduction of 
the mean dipole size, making the medium more transpar- 
ent. Correspondingly, we should expect that the analysis 
of ALICE data, performed in [7[, should have underes- 
timated the medium density, i.e. the parameter qq. We 
will return to this problem in Sect. IIV Al 

Unfortunately, Eq. (|24|) has an analytic solution only 
if e(l) is constant, which is not the case here. Then one 
should solve the equation numerically, so the simplicity 
of such a heuristic model is lost. In these circumstances 
it is worth switching to the rigorous quantum-mechanical 
description of the evolution and attenuation of a dipole 
in an absorptive medium, and employ the path-integral 
approach [421 ]. 



B. Path-integral technique 

Evolution and attenuation of a qq dipole propagating 
through an absorptive medium, starting from the trans- 
verse qq separation r\ at a point with longitudinal co- 
ordinate 1 1 and evolving its size up to r*2 at the point 
I2, is given by a sum over all possible trajectories of q 
and q. The resulting survival probability amplitude has 
the form of a light-cone Green function Gq q (h, r*i; I2, fa), 
which satisfies the two-dimensional Schrodinger equation 

SUM, 



1 Factor 2 was missed in 0] 



d/2 2pT a (1 — a) 

= id(l 2 - h) S(r 2 - n) 



Vq q (h,r 2 ) 



Gq q (h,f 'i;h,r 2 ) 
(25) 
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and the boundary conditions, 

Gg q (ll,fi;l 2 ,f2) 



The second term in square brackets in (|2"5l plays role 
of the kinetic energy in Schrodinger equation, while the 
imaginary part of the light-cone potential Vq q (l2,r2) in 
(1251) is responsible for absorption in the medium. Ac- 
cording to Eqs. CEHI-CG 



The real part of the potential describes the nonpertur- 
bative interaction between q and q in the dipole [22], |42j . 
= S(r 2 — ri)\ It is questionable, however, whether such a binding po- 

tential should be considered within a hot, maybe de- 
= o (26) confined medium. Therefore we will treat the qq as free 

noninteracting partons, like in the previous Sect. IIII Al 
The real potential should not affect much the dipole evo- 
lution on the initial perturbative stage of development. 



ImV q : g (l,r) = -~q(l) 



In the case of a constant medium density q(l) = q = 
(27) const, Eq. ([25]) allows an analytic solution |42J, 



Gqq(h,ri;l 2 ,r2) 



2-ki sin(o;AZ) 



exp • 



i ^ 



2sin(wAZ) 



{r\ + r\) cos(wAZ) - 2r*i • f 2 



exp 



2pTCt(l — a) 



(28) 



where Al = l 2 — h and 
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2 ptol{1 — - a) ' 

i 

--p T a(l - a)q. 



(29) 



It was demonstrated in [42J, |43| that in the general case 
of arbitrarily varying medium density the Green function 
still retains the oscillatory form, 



Gq q {h,r \\l 2 ,r 2 ) = fi(li,Z 2 )exp 



im 2 q M 



(30) 



2pTCt(l — a) 

x exp {Tx&./a) r\ + T 2 (h,h) r\ + T 3 (h,h) ?i • r 2 } . 

All information about the dipole absorption rate varying 
along its trajectory, q(l), is contained in the coefficients 



£1 and Tj. They are calculated numerically by slicing the 
medium into layers, which are sufficiently thin to keep 
q, given by (|20|) . constant within each of them. Then 
one can calculate the coefficients 0(Z X , Z 2 ) and Ti(li,l 2 ), 
employing the recurrence relations between subsequent 
layers, derived in [42l l43j. 

As in the previous section, we rely on the Berger model 
i40j, which assumes equal sharing of the pion light-cone 
momentum between q and q, i.e. a = 1/2. Since the 
Green function is projected into the pion light-cone wave 
function, we have to keep a — 1/2 for the dipole as well. 

Now we are in a position to write a rigorous quantum- 
mechanical extension of the simplified model (|23|) . for the 
suppression factor Rab (b) of hadrons produced with high 
Pt in a hard process in a collision of nuclei A and B with 
relative impact parameter b, 



J 



2^ 



Rab (b, pt) 



fd*TT A (T)T B (b-T)f% 




J da J d 2 rid 2 r2^l(r2,a)Gq q (h,r 1 -,l2,r2)^ i n{ri,a) 
o 



Tab (b) 



J da J d 2 r^l(f 2 :a)^in(ri,a) 



T AB (b) |* fc (0)| 



, d 2 TT A (r)T B (b-T) 



drr* h (r) Gq q (0,0; 



(31) 



where Tab — J d 2 rTA(b)TB(b — r); </> is the azimuthal relative to the impact vector b of the collision. l max is 
angle of the dipole trajectory in impact parameter plane, any distance, which should be much longer that the ex- 
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tent of the medium. Its particular length is unimportant, 
because the Green function in vacuum is just a phase. 
Notice that all information about the dipole trajectory, 
including the </>-, f- and 6-dependences, is contained in 
the Green function. In (|3"Tj) we fixed r = 0, because ac- 
cording to the solution Eq. (|22|) , it does make a difference 
whether tq is as small as 1/kx or zero. It is also worth 
noting that we start the evolution of the Green func- 
tion Gq q {l\ 1 fx] I2, ^2) at l\ — 0, i.e. at the point of the 
hard collision, but assume that the absorptive imaginary 
part of the potential Eq. (1231) is zero at I < Iq (compare 
with (j2"0")) '). Although the dipole does not exist at I < l p , 
the parton virtuality is steeply falling, governed by the 
same equation (|2ip . and the dipole is produced at I = l p 
with an enlarged separation, the same as if it had started 
evolution at t = 0. We remind that the mean produc- 
tion length (lp) is different for quark and gluon jets as is 
demonstrated in Fig. [SJ Therefore the numerator in (|31[) 
is calculated separately for quark and gluon jets and then 
summed with the weights given by Eq. (fTST) . 

During the short path from I = Iq to I = l p (if l p > 
Iq) the parton experiences multiple interactions, which 
induce extra radiation of gluons and additional loss of 
energy [27j|, 



i P 1 

AE = 9(Z P - Jo) J dl J dl' q(l'). 



(32) 



that CT is the source of the rising px dependence of 
Raa observed in the ALICE experiment 41 1 . 

Encouraged by the success of the simple model, we 
perform here full calculations employing the path integral 
method, Eq. (13T1) . The results for central (0 — 5%) lead- 
lead collisions at yfs = 2.76 TeV are shown by the dashed 
curve in Fig. [51 compared with new data from the ALICE 
[HI and CMS [4{| experiments, extended to higher 
values of px than those in [4l[. The only free parameter, 



Pb-Pb Vs NN = 2.76 TeV 
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Although this is a small correction, we included it in 
the calculations by making a proper shift of the variable 
Zh in the fragmentation function. 



IV. COMPARISON WITH DATA 

Numerous results of new precise measurements at 
RHIC and LHC have been released recently. They allow 
to perform stringent tests of the contemporary models of 
in-medium hadronization. 



A. Quenching of high-pr hadrons 

The comparison of RAA(b, px) calculated within the 
simple model described in Sect. IIII Al with data from 
the ALICE experiment [4l[ at *fs = 2.76 TeV, was per- 
formed in 0. While the absolute value of Raa is ad- 
justable, its rising py-dependence originates from the re- 
duction of the mean dipole size with px, in accordance 
with Eq. (l22l) . and due to Lorentz dilation of the dipole 
size expansion. As a result, the medium becomes more 
transparent for more energetic dipoles in accordance with 
the effect of color transparency (CT), which makes the 
medium more transparent for smaller dipoles. An anal- 
ogous rising energy dependence of medium transparency 
was predicted and observed for virtual photoproduction 
of vector mesons on nuclei 1311 . So it was concluded in 



FIG. 8: (Color online) The suppression factor Raa for cen- 
tral (0 - 5%) lead-lead collisions at y/s = 2.76 TeV. The 
dashed line is calculated within the path-integral approach, 
Eq. (|31|) with the space- and time dependent transport coeffi- 
cient Eq. ()20[) . where the adjusted parameter qo — 2 GeV 2 / mi. 
The solid curve also includes the effects of initial state interac- 
tions in nuclear collisions 6, 51], as is described in Scct. HVB"] 
Data for Raa are from the ALICE [3 and CMS gaUIl 
experiments. 

the maximal transport coefficient defined in Eq. (|20p , was 
adjusted to the data and fixed at % = 2 GeV 2 /fm for all 
further calculations for lead-lead collisions at this energy. 

Notice that this value of qo is about twice larger than 
what was found in 0] within the simple model described 
above in Sect. IIII Al As we commented, this should be 
expected, since Eq. (|2"Tj) is lacking the color filtering ef- 
fect added in (j2"4"|) or inherited from the Green function 
equation (j2"o| . 

While our calculations describe well the data at high 
Pt > 6 GeV, the region of smaller px is apparently dom- 
inated by a thermal mechanisms of hadron production, 
as is confirmed by the large elliptic flow observed in this 
region (see below). This is why we do not try to explain 
data in this region by the perturbative dynamics. 

The variation of the suppression factor Raa(pt, b) 
with impact parameter of collision was also calculated 
with the Eq. (|3"Tj) . The results plotted by dashed curves 
are compared with data taken at different centralities of 
collision by the ALICE experiment [48| in Fig. |H1 and by 
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the CMS experiment [H, [Ho| in Fig. [TDJ In all cases we 
observe good agreement. 



Towards large xt- new constraints from 
energy conservation 
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FIG. 9: (Color online) Centrality dependence of the suppres- 
sion factor Raa{pt ,b) measured in the ALICE experiment 
pH ]. The intervals of centrality are indicated in the plot. The 
meaning of the curves is the same as in Fig. [8] 
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FIG. 10: (Color online) The same as in Fig. [9] with data from 
the CMS experiment g! Hi- 



lt was stressed in [H [oT| that energy conservation may 
become an issue upon approaching the kinematic bound 
of either large Feynman xl = Xf = %Pl/\/~s, or trans- 
verse fractional momentum xt = 2px/\/s. Apparently, 
in such a kinematic domain any initial state interac- 
tions (ISI) leading to energy dissipation, should result 
in a suppressed production rate ofparticles with large 
xl,t- Indeed, as it was stressed in [6(, every reaction ex- 
perimentally studied so far, with any particle (hadrons, 
Drell-Yan dileptons, charmonium, etc.) produced with 
large xl, exhibits nuclear suppression. Observation of 
such a suppression in high xt processes is more difficult, 
because the cross sections steeply fall with px, and avail- 
able statistics may not be sufficient (see however below 
Fig. [TTj) . As in [51j , we apply to high-p-r production ex- 
actly the same model developed in (y, |52| for large Xl , 
and with the same parameters. 

Since initial state multiple collisions suppress the pro- 
duction rate of leading particles, we assume that every 
collision brings in a suppression factor U(£) [6[, where 
£ = \J x\ + x T This factor [/(£) should cause a strong 
(for heavy nuclei) suppression at £ — > 1, although at the 
same time some enhancement at small £ due to the feed 
down from higher £. This is because energy conservation 
does not lead to disappearance (absorption) of particles, 
but only to their re-distribution in £. 

Since at £ — > 1 the kinematics of an inelastic collision 
corresponds to no particle production within the rapidity 
interval Ay ~ — ln(l— £), the suppression factor U (£) can 
also be treated as survival probability of a large rapidity 
gap, which is Sudakov suppressed. The mean number 
(rig (Ay)) of gluons radiated in the rapidity interval Ay is 
related to the height of the plateau in the gluon spectrum, 
(rig (Ay)) = Ay drigj dy. Then, the Sudakov factor reads, 



U(0 = (1 - £) dr ' 



,/ d v 



(33) 



The height of the gluon plateau was estimated in [53[ as, 



drig 
dy 



3a 



iir 



hi , . , 

n \ A QCD 



(34) 



The value of a s was fitted in [53[ at a s = 0.45, using data 
on pion multiplicity in e + e~ annihilation. This leads 
with a good accuracy to dn g /dy s» 1, i.e. the Sudakov 
factor, 



u(o = i - e 



(35) 



Although QCD factorization is expected to be broken 
by ISI in pA collisions at large £, we will rely on the effec- 
tive factorization formula, Eq. (|15[) , where we replace the 
proton parton distribution function by a nuclear modi- 
fied one F i / p (x il Q 2 ) => F^(x t , Q 2 , b). In the case of 
nuclear collisions we do this modification for the bound 



12 



nucleons in both nuclei. Relying on the suppression fac 



tor Eq. (|35l) . and applying the AGK cutting rules [54 1 
with the Glauber weight factors, one arrives at the ISI 
modified parton distribution function of the proton in a 
pA collision at impact parameter &, 



^i/p (^* 



Q\b) = CF l/p { Xl ,Q 2 ) 



fT A (b)l 



(1-6 1 



fT A (b)] 



(36) 



Here a e f f is the effective hadronic cross section control- 
ling multiple interactions. It is reduced by Gribov inelas- 
tic shadowing, which makes the nuclear medium much 
more transparent. The effective cross section was evalu- 
ated in 0, [55| at about a e f / « 20 mb. The normalization 
factor C in Eq. ([31)]) is fixed by the Gottfried sum rule, be- 
cause the number of valence quarks, dominating at large 
£, should be unchanged. 

With the parton distribution functions Eq. (|3"6"]) mod- 
ified by ISI one achieves a good parameter-free descrip- 
tion of available data at large xl [1, [H|. Moreover, these 
corrections may be important at large pt, in particu- 
lar in the RHIC energy range, where xt reaches val- 
ues of 0.2 — 0.3. Notice that the real values of xt, es- 
sential for energy conservation, are significantly higher, 
xt = xr/zfi, so it reaches values of 0.3 — 0.5 at RHIC, 
and about 0.4 at LHC (p T = 200 GeV). 

Comparison with PHENIX data [56| for neutral pion 
production in central d-Au collisions is presented in 
Fig. 111! There is a firm indication in data that the pre- 



In nuclear collisions the ISI effects are calculated sim- 
ilarly, using the modified parton distribution functions 
Eq. ([55)1 for nucleons in both colliding nuclei. The re- 
sulting complementary suppression reduces Raa(pt) at 
large xt- This is demonstrated in Fig. [12] If we just re- 
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FIG. 12: (Color online) Nuclear attenuation factor Raa(pt) 
for neutral pions produced in central gold-gold collisions at 
x/s = 200 GeV. The solid and dashed line are calculated 
with or without ISI corrections. PHENIX data are from [57l ] 
(triangles) and [H3] (squares). 
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peated the same calculations as done above for the LHC, 
we would get Raa steeply rising with px, as is depicted 
by the dashed curve. However the ISI effects and energy 
conservation impose a sizable additional suppression, as 
is shown by the solid curve. Apparently this improves 
the agreement with data. Notice that the only fitted pa- 
rameter, the transport coefficient, should be re-adjusted 
and was found to be go = 1.6 GeV 2 / fm at this energy, 
smaller than at y/s — 2.76 TeV, as expected. 

Keeping this parameter fixed we can calculate other 
observables for gold-gold collisions at y 7 ! = 200 GeV. 
Fig. H31 presents our results for the suppression of neutral 
pions at different collision centralities, in comparison 
with PHENIX data (53 . 



FIG. 11: (Color online) Nuclear attenuation factor RdAu(pT) 
as function of pr of tv° mesons produced in central (0 — 20%) 
d-Au collisions at ^/s = 200 GeV and r\ = 0. The solid and 
dashed curves represent the model predictions calculated with 
and without the ISI corrections. Isotopic effect is included. 
The data are from the PHENIX experiment [5a ]. 

dieted strong nuclear suppression at large xt is indeed 
observed. Still, the statistical evidence of the effect needs 
to be enhanced. 



One can access even larger values of xt = %Pt / y/s, 
by either increasing px or going down in energy. In 
both cases the effects of ISI should be stronger. Indeed, 
data for Raa(pt) in gold-gold collisions at = 62 GeV 
plotted in Fig. [TJ] show a falling, rather than rising 
PT-dependencc. Our calculations using Eq. ([3Tj) again 
demonstrate good agreement. Of course, the hot medium 
properties changed, so we had to re-adjust the parameter 
go = 1.2GeV 2 /fm. 
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FIG. 13: (Color online) Centrality dependence of the suppres- 
sion factor Raa^Pt,V) measured in the PHENIX experiment 
[53| in gold-gold collisions at yfs = 200 GeV. The intervals of 
centrality are indicated in the plot. 



C. Azimuthal anisotropy 

Experimental observation of a suppression of high-p^ 
hadrons escaping from the dense medium means that not 
the whole volume of the medium contributes. Indepen- 
dently of the suppression mechanism (except for the con- 
tribution of ISI) , this means that the longer is the path- 
length of propagation in the medium (in two-dimensional 
transverse plane), the stronger is the suppression. Thus, 
one can conclude that the direction of propagation nor- 
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FIG. 15: (Color online) Nuclear factor Raa(pt) for charge 
hadrons in lead-lead collisions at \/s = 2.76 TeV at differ- 
ent centralities. ALICE data {6lj ] and our calculations with 
Eq. (|3ip are divided into two classes: In plane (37r/4 < \<f>\ < n 
and |0| < 7r/4), and Out-of-plane (tt/4 < \<f>\ < 37r/4). 
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FIG. 14: (Color online) The same as in Fig. \T3\ but at y/s ■ 
62 GeV. Data are from \M- 



mally to the medium surface is preferable. For a non- 
central collision with an almond shape of the intersection 
area this should lead to an azimuthal asymmetry of high- 
Pt hadron production. Although this conclusion should 
be valid for any mechanism of suppression, the magni- 
tude of the asymmetry is of course model dependent. 

Our main result for high-py hadron suppression, 
Eq. (|3"Tj) . can be tested by comparison of the predicted 
azimuthal angle dependence with data. Following data 
(6l| , we split the integration over <j> in d3T]) into two inter- 
vals: (i) |0| > 37r/4plus 101 < tt/4; (ii) tt/4 < \<j>\ < 3tt/4. 
These two contributions are called in- and out-of-plane 
respectively. As expected, in-plane events are less sup- 
pressed compared with out-of-plane, as is depicted by the 
upper and bottom curves in Fig. [15j Our results agree 
well with ALICE data 61] at pt > 6 GeV and at all 
measured centralities. 

Similar measurements of Raa were performed in the 
PHENIX experiment for the In- and Out-of the scattering 
plane events 57|. The results for centrality 20-30% are 
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plotted in Fig. [T5] in comparison with our calculations. 
Notice that both Figs . [Lol and [TBI show that the transition 
from the hydrodynamic to perturbative regimes occur 
for In- plane events with a delay, at higher px- This is 
natural, since the hydrodynamic flow is much stronger, 
correspondingly the cross section is larger. 



FIG. 16: (Color online) Raa(pt) for charge hadrons in gold- 
gold collisions at y/s = 200 GeV at centralities — 20%. 
PHENIX data (57J and our calculations with Eq. ([31) are di- 
vided into two classes of events: In-plane (ll7r/12 < \cj>\ < n 
and |0| < 7T/12), and Out-of-plane (5tt/12 < \<j>\ < 7tt/12). 



Usually data on azimuthal asymmetry of particle pro- 
duction are presented in terms of the second moment of 
the (^-distribution, v 2 = (cos(20)). We can calculate it 
with a slight modification of Eq. (j3"Tj) , 



/ d 2 TT A (T)T B (b~ f) J d0cos(20) 



v 2 (pT,b) 



l maxi 1 J 





JdHT A (r)T B (b-r)J 
o 



/ drrV h {r)Gq q (0,0;l max ,r) 
o 



(37) 



We remind that the Green function implicitly depends on 
the impact parameters b and r, and on the trajectory of 
the dipole in the hot matter, and that these dependences 
are contained in the functions Q(li,l2) and Y^Zi,^) in 
Eq. ([30]). 

In order to compare with data on i>2 (pt), one should 
integrate the numerator and denominator in (|37|) over 
the intervals of impact parameter from b m i n to b max , 
which correspond to the measured intervals of central- 
ity. Our results are compared with ALICE data [6l| and 
with CMS data in Figs. Q2] and QH respectively. In 
all cases we observe good agreement. 

Naturally, as with data on Raa, our pQCD calcu- 
lations for V2(pt) grossly underestimate data at small 
Pt ^ 6 GeV. Remarkably, the transition to the pQCD 
regime occurs for V2(pr) at the same pt, as for Raa{pt)- 
This confirms the presence of two different mechanisms: 
the dominant hydrodynamic mechanism of elliptic flow, 
providing a large and rising with px anisotropy «2(Pr)> 
which abruptly switches to the regime of pQCD, having a 
much smaller azimuthal anisotropy. So it is not acciden- 
tal that both Raa(pt) and i>2 (pt) swiftly change their 
behavior at the same value of px- 

We also compare our results for the azimuthal 
anisotropy with RHIC data. Although we included the 
effects of ISI suppressing Raa at large px, these correc- 
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FIG. 17: (Color online) ALICE data [H for azimuthal 
anisotropy, V2, vs pr for charge hadron production in lead- 
lead collisions at mid rapidity, at y/s = 2.76 TeV and at dif- 
ferent centralities indicated in the figure. The curves present 
the results of calculation with Eq. ()37l) . 
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FIG. 18: (Color online) The same as in Fig. 1171 but displaying 
data from the CMS experiment [oil ]. 



tions mostly cancel in the ratio, Eq. (|57|) . Our calcu- 
lations agree well with PHENIX data for the azimuthal 
asymmetry for 7r° production at s/s = 200 GeV and at 
mid rapidity, as is demonstrated in Fig. 1191 
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FIG. 19: (Color online) The same as in Fig. 1171 but displaying 
data from the PHENIX experiment [6(| . 



V. SUMMARY AND PROSPECTIVES 

This paper attempts at a quantitative understanding 
of experimentally observed strong attenuation of hadrons 
inclusively produced with large transverse momenta in 
heavy ion collisions, based on contemporary models for 
space-time development of hadronization. The improved 
quality of data from RHIC and new high-statistics data 
from LHC makes the test of models more challenging 
and also more decisive. In particular, the popular energy 
loss scenario, based on the unjustified assumption of long 
production length, experiences difficulties explaining the 
new data. It also fails to explain the data from the HER- 
MES experiment at HERA for leading hadron production 
in semi-inclusive DIS, which provides a sensitive testing 
ground for in-medium hadronization models. 

Here we presented an alternative mechanism of sup- 
pression of high-pT" hadrons inclusively produced in 
heavy ion collisions. First of all, we investigated the 
space-time development of gluon radiation and energy 
dissipation by a highly virtual parton produced in a high- 
er process. The key point of this consideration is an 
uniquely high initial virtuality of such a parton, which is 
of the order of its energy. Thus, increasing the jet energy 
one simultaneously enhances the hardness of the process, 
intensifying the dissipation of the energy. For this rea- 
son energy conservation imposes tough constraints on the 
production length of leading hadrons, which does not rise 
with pt, but remains short, as is demonstrated in Fig. [5] 

It is worth emphasizing that the shortness of the pro- 
duction length is a specific feature of inclusive hadron 
production. However, formation of a high-p^ jet is char- 
acterized by a much longer time scale. Indeed, the convo- 
lution Eq. (TT5]) of steeply falling cross section of parton- 
parton scattering with parton distribution and fragmen- 
tation functions strongly enhances the contribution of 
large Zh, as is demonstrated in Fig. [7] In this case energy 
conservation is an important issue and restricts the pro- 
duction length. As for fully reconstructed high-pf jets, 
they are characterized by a very different space-time de- 
velopment, and the medium-induced energy loss may be 
indeed an important source of the observed jet suppres- 
sion. 

With the production of a colorless hadronic state 
(dipole) the gluon radiation and energy loss cease, and 
the main reason for suppression of the production rate 
becomes the survival probability of the pre-hadron prop- 
agating through the dense matter. Apparently, at larger 
Pt the expansion of the initially small initial size of the 
dipole slows down due to Lorentz time dilation. Here the 
effect of color transparency is at work, and the medium 
becomes more transparent for more energetic dipoles. A 
rising nuclear suppression factor Raa{pt) is indeed ob- 
served at LHC, in good accord with the effect of color 
transparency. 

In this paper we performed calculations based on the 
most strict quantum-mechanical description of space- 
time development and attenuation of color-dipoles prop- 
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agating through a medium, the path-integral method. 
We calculated and compared with data the suppression 
factor RaaQ>) for different ccntralities and energies of 
collision, ranging from y/s = 62 GeV to 2.67 TeV. Ad- 
ditional suppression arising from initial state multiple 
interactions, important at large xt or x^, was also in- 
cluded. The related corrections are found to be impor- 
tant at xt > 0.1, where they slow down the rise of Raa 
and even turn it into a falling p^-dependence. We also 
calculated the azimuthal anisotropy of hadron produc- 
tion, which reflects the asymmetric shape of the overlap 
area of the colliding nuclei. 

In all cases we reached good agreement with available 
data from RHIC and LHC at high pr- The only adjusted 
parameter, the maximal value of the transport coefficient, 
Eq. flU}, was found to be q Q = 2 GeV 2 /fm, 1.6GeV 2 /fm 
and 1.2 GeV 2 /fm at y/s = 2.76 TeV, 200 GeV and 62 GeV 
respectively, for heavy nuclei, lead and gold. 

Our results reproduce quite well data at pt > 6 GeV. 
However, the observed Raa(j>t) and ^(pt) expose quite 
a different behavior towards smaller px, steeply rising 
and shaping a bump. We relate this to an interplay of 
two mechanisms of hadron production: (i) evaporation of 
hadrons from the created hot medium controlled by hy- 
drodynamics; (ii) perturbative QCD mechanism of high- 



Pt production of hadrons, which propagate and attenu- 
ate in the hot medium. The abrupt transition between 
the two mechanisms causes distinct minima in Raa{pt) 
and in V2(pt), both observed at the same values of pt- 
We plan to work on combining the two mechanisms, aim- 
ing to describe data in the full range oipT- 
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